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Abstract 

This  paper  presents  a  flight  controller  for  a  tail¬ 
less  aircraft  with  a  large  suite  of  conventional  and 
unconventional  control  effectors.  The  controller 
structure  is  modular  to  take  advantages  of  individual 
technologies  from  the  areas  of  plant  parameter  esti¬ 
mation,  control  allocation,  and  robust  feedback  con¬ 
trol.  Linear  models  generated  off-line  provide  plant 
parameter  estimates  for  control.  Dynamic  inversion 
control  provides  direct  satisfaction  of  flying  qualities 
requirements  in  the  presence  of  uncertainties.  The 
focus  of  this  paper,  however,  is  control  allocation. 
Control  allocation  is  posed  as  constrained  param¬ 
eter  optimization  to  minimize  an  objective  that  is 
a  function  of  the  control  surface  deflections.  The 
control  law  is  decomposed  into  a  sequence  of  prior¬ 
itized  partitions,  and  additional  optimization  vari¬ 
ables  scale  the  control  partitions  to  to  provide  opti¬ 
mal  command  limiting  which  prevent  actuator  sa^ 
uration.  Analysis  shows  that  appropriate  prioriti¬ 
zation  of  dynamic  inversion  control  laws  provides 
graceful  command  and  loop  response  degradation  for 
unachievable  commands.  Preliminary  simulation  re¬ 
sults  show  that  command  variable  response  remains 
decoupled  for  unachievable  commands  whUe  other 
command  limiting  methods  may  result  in  unaccept¬ 
able  coupled  response. 

1  Introduction 

The  next  generation  United  States  Air  Force  fighter 
aircraft  will  more  than  likely  have  a  reduced  verti¬ 
cal  tail  or  no  vertical  tail  at  all.  Radar  cross  sec¬ 
tion  and  structural  weight  reduction  benefits  have 
influenced  these  tailless  aircraft  configurations  [1]. 
Tailless  configurations,  however,  present  an  inter¬ 
esting  challenge  from  a  stability  and  control  per¬ 
spective  [2].  Vertical  tails  provide  directional  sta¬ 
bility  and  contain  the  rudder  which  is  typically  the 
primary  directional  control  effector.  Lack  of  a  vet- 
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tical  tail  undoubtedly  relaxes  directional  stability 
and  reduces  directional  control  power.  Advanced 
fighter  configurations  may  also  have  relaxed  longitu¬ 
dinal  stability  for  reduced  supersonic  trim  drag  and 
enhanced  maneuverability  which  further  challenges 
the  flight  control  system.  Unconventional  control  ef¬ 
fectors  that  provide  multi-axis  forces  and  moments 
may  be  added  to  overcome  control  power  deficiencies 
[3,  1].  The  unconventional  control  effectors  must  be 
allocated  with  other  conventional  effectors  to  main¬ 
tain  performance  of  current  fighters.  Flight  con¬ 
trol  algorithms  must  augment  the  releixed  stability 
of  these  configurations  and  manage  the  distribution 
of  many  control  effectors.  Flight  control  technology 
must  therefore  pave  the  way  for  successful  tailless 
aircraft  programs. 

A  modular  control  structure  for  flight  control  of  a 
tailless  aircraft  configuration  is  developed  to  exploit 
technological  advances  in  feedback  control,  param¬ 
eter  estimation  and  control  allocation.  The  flight 
controller  consists  of  feedback  and  feedforward  aug¬ 
mentation,  control  allocation,  and  parameter  esti¬ 
mation  modules.  The  parameter  estimation  module 
provides  stability  and  control  derivative  estimates 
to  the  control  allocation  and  augmentation  modules. 
The  augmentation  module  provides  flight  condition 
adaptation,' desired  flying  qualities,  and  robustness 
to  uncertainties.  The  control  allocation  module  allo¬ 
cates  the  redundant  limited  control  effectors  to  opti¬ 
mize  an  objective,  achieve  maximum  control  power, 
and  limit  pilot  commands  to  reduce  the  risk  of  de¬ 
parture. 

The  modular  flight  control  structure  has  been 
initially  applied  to  a  tailless  fighter  configuration 
[3].  The  flight  controller  included  an  inner  equal¬ 
ization/outer  robust  loop  augmentation  module,  an 
off-line  generated  database  plant  estimation  module, 
and  simple  ganging  control  allocation  module  that 
provided  no  command  limiting.  The  flight  control 
design  was  successfully  demonstrated  for  flight  of  a 
tailless  aircraft  configuration  with  minimal  uncon- 
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ventional  control  effectors  in  a  limited  envelope. 

This  paper  presents  extensions  of  past  work  [3] 
by  applying  the  modular  control  structure  to  flight 
control  of  a  tailless  fighter  with  a  larger  effector  suite 
than  previous  work.  The  most  significant  extension 
to  past  work,  however,  is  to  explore  the  full  con¬ 
trol  allocation  capability  of  the  tailless  fighter  air¬ 
craft.  Optimal  objective  control  allocation  and  pri¬ 
oritized  control  command  limiting  are  formulated  as 
a  constrained  parameter  optimization  problem.  An 
objective  function  of  the  control  deflections  is  con¬ 
structed,  and  achievable  commands  are  realized  with 
control  deflections  that  minimize  the  objective.  Un¬ 
achievable  commands  are  limited  and  unstable  axes 
prioritized  by  decomposing  the  control  law  into  a 
prioritized  sequence  of  commands.  Scalar  multipli¬ 
ers  from  the  parameter  optimization  solution  indi¬ 
cate  whether  the  command  is  achievable  or  not  and 
provide  scalings  for  the  sequence  of  commands  that 
provide  optimal  command  limiting  and  axes  priori¬ 
tization. 

2  Model  of  Aircraft  with 
Redundant  Effectors 

The  aircraft  nonlinear  equations  of  motion  are  well 
documented  [4,  5]  and  can  be  represented  by  a  sys¬ 
tem  of  ordinary  differential  equations  of  the  follow¬ 
ing  form 

i(0  =  /(x(0.«(0)  (1) 

where  x  6  IR”*  is  the  state  and  ti  G  IR”*  is  the 
control.  At  some  arbitrary  time  along  an  arbitrary 
trajectory,  the  aircraft  dynamics  are  approximated 
by  the  following  set  of  time- varying  differential  equa¬ 
tions 

i(0  =  .4(0x(0  +  B(0«(0  +  6(0  (2) 

where  A{t)  are  the  stability  derivatives,  B(t)  are  the 
control  derivatives,  and  b(t)  is  a  bias  that  includes 
higher  order  terms  and  accelerations  due  to  non¬ 
equilibrium  or  off-trim  conditions.  These  equations 
are  nonlinear  in  time  and  non-homogeneous  due  to 
the  bias. 

A  set  of  commanded  variables  y  E  IR”^  is  defined 
by  the  user  as  a  linear  function  of  the  state 

y(0  =  (3) 

rank(Cy(^))  =  riy  . 

Note  that  eq.(3)  is  nonlinear  in  time.  As  shown  in 
[6],  the  following  non-unique  state  transformation 
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which  transforms  the  system  defined  by  eqs.(2)  and 
(3)  into  the  following  convenient  form 


r  m  1 

_  r  ■Azzii) 

Azy{t)  1 

[  <i)  1 

1  m  \ 

■  [AzW 

-^yy(0  J 

1  y(0  J 

.  [5,(0 
+  [  By{t) 

]«(<)+! 

[  m ' 

L  ^y(^)  . 

(5) 


Also  shown  in  [6],  the  control  dimension  may  be  re¬ 
duced  for  design  by  introducing  a  control  allocation 
function  p  if  the  system  in  eq,(5)  is  redundant 
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where  control  redundancy  is  defined  by  the  following 
conditions 

rank{by,)  =  lVt  =  l,-;,m  (7) 

rank{By)  =  Uy 

m  >  Uy 

By  =  [  by^  •••  by^  ]  . 

\ 

The  square  system  in  eq-(6)  provides  the  design 
model  for  the  flight  controller  given  in  the  follow¬ 
ing  section. 


3  Controller  Structure 

The  proposed  flight  controller  for  the  aircraft  model 
in  section  2  has  a  modular  structure  that  benefits 
from  independent  technology  advances  in  feedback 
control,  control  allocation  and  plant  parameter  es¬ 
timation.  The  controller  structure  is  given  in  Fig. 
1.  The  Feedback  and  Feedforward  Augmentation 
(FFA)  module  inverts  the  output  dynamics  and  com¬ 
mands  appropriate  dynamics  to  provide  desired  sta¬ 
bility  and  flying  qualities  properties  in  a  manner 
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that  is  robust  to  uncertainties.  The  Control  Allo¬ 
cation  (CA)  module  provides  actuator  commands 
that  optimize  an  objective,  utilize  maximum  control 
authority,  and  limit  unachievable  commands.  The 
plant  MODEL  module  provides  stability  and  con¬ 
trol  derivative  estimates  of  the  AirCraft  (A/C)  to 
the  FFA  and  CA  modules. 

3.1  Plant  Model 

Currently,  off-line  generated  linear  models  are  used 
to  supply  aircraft  stability  and  control  derivatives 
to  the  augmentation  and  control  allocation  modules. 
The  nonlinear  aircraft  model  moments  are  balanced 
for  wings-level  flight  at  a  given  Mach  number,  al¬ 
titude,  and  angle-of-attack  to  establish  a  nominal 
flight  condition.  Perturbation al  techniques  are  then 
used  to  generate  a  linear  model  of  the  aircraft  dy¬ 
namics  near  the  nominal  flight  condition  of  the  fol¬ 
lowing  form  which  is  identical  to  the  model  in  eq.(2) 
at  some  fixed  time 

x(<)  =  Ax{t)  -1-  Bu{t)  +  b.  (8) 

Note  the  bias  term  b  that  results  due  to  untrimmed 
forces. 

3.2  Feedback  and  Feedforward 
Augmentation 

Dyneunic  inversion  and  proportional/integral  com¬ 
pensation  is  used  for  robust  feedback  stability  aug¬ 
mentation,  and  feedforward  compensation  is  used 
for  command  augmentation  [7,  6].  Dynamic  in¬ 
version  inverts  the  output  dynamics  and  propor¬ 
tional/integral  feedback  generates  a  desired  robust 
loop  shape.  Proportional  gains  are  used  in  the  feed¬ 
forward  path  to  obtain  desirable  flying  qualities. 

Derivation  of  the  control  law  begins  with  the 
plant  model  in  eq.(6).  Recall  that  this  model  is 
square  since  the  number  of  controls  dy  equals  the 
number  of  outputs  y.  The  dynamic  inversion  con¬ 
trol  law  is  given  by  [6] 

dy**  =  ~AyzZ{t)  —  Ayyy(t')  —  Jy  t>  •  (9) 

Additional  control  is  accomplished  through  input 
V.  The  following  proportional-plus-integral  control 
with  feedforward  elements  around  the  integrators 
has  been  shown  [7,  6]  to  provide  satisfactory  robust¬ 
ness  and  flying  qualities  for  fighter  aircraft 

V  =  Wclj  —  U>cy  -b  We/ey**""*  (10) 

Xi  =  -w*/<y-l-a;*/,y*’”'*. 


The  feedforward  gain  is  fc,  and  the  proportional  and 
integral  gains  are  functions  of  cjc  and  /*.  The  closed 
loop  system  is  given  by 
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A(z.y,r.-,y*"“')  =  5.p(dy"**).  (11) 


Note  that  if  is  sufficiently  small  [6],  then  the  sys¬ 
tem  in  cq.(ll)  is  approximately  linear.  Further  for 
Bz  =  0,  exarhination  of  the  characteristic  equation 
zeros  may  be  considered  for  nominal  stability  anal¬ 
ysis 


det(AJn.+2n,-^l)=  (12) 

det  (A/n.  ~  (A"  +  o/cA  +  =  0 . 

So  if  det(A/n,  -  A^z)  has  no  zeros  in  the  right  half  of 
the  complex  plane,  nominal  stability  is  insured  with 
a  proper  choice  of  u)c  and  /*. 

The  command/output  channels  are  decoupled 
due  to  the  noninteracting  property  of  the  dynamic 
inversion  control,  and  the  individual  channel  trans¬ 
fer  functions  are  given  by  the  following 

Vi  _  '^tfcS->rUilSi  /jgx 

where  j  =  The  command  channel  re¬ 

sponse  is  therefore  directly  related  to  the  control 
parameters.  So  nominal  stability  and  performance 
are  readily  satisfied  with  a  proportional-plus-integral 
control  structure.  However,  the  proportional  and 
integral  feedback  gains  must  also  account  for  sys¬ 
tem  uncertainties.  Loopshaping  [6]  provides  a  simple 
procedure  for  choosing  robust  feedback  gains  since 
the  noninteracting  property  of  the  dynamic  inversion 
decouples  the  MIMO  aircraft. 


3,3  Control  Allocation 

It  is  assumed  that  the  control  law  in  eqs.(9)  and 
(10)  provides  desired  loop  and  command  response 
properties,  however  the  reduced  dimension  controls 
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must  be  resolved  into  actual  control  effector 
commands.  The  control  allocation  function  in  eq.(6) 
provides  the  mapping  from  reduced  dimension  con- 
trols  to  actual  control  commands  u  = 

The  control  allocation  problem  has  been  stated 
many  times  in  various  forms  [8,  6,  9].  The  unlimited 
control  allocation  problem  originates  from  the  con¬ 
trol  dimension  reduction  in  eq.(6).  The  unlimited 
control  allocation  problem,  expressed  analytically, 
is  to  find  a  function  p  such  that  Byp{dy^*)  =  dy^* 
for  all  dy^*.  For  redundant  control  effector  suites, 
there  are  many  functions  that  solve  the  control  allo¬ 
cation  problem.  The  effector  redundancy  maybe  ex¬ 
ploited  to  satisfy  additional  objectives  such  as  min¬ 
imal  radar  cross  section,  drag,  or  wing  loads.  How¬ 
ever,  since  there  are  limits  on  the  control  effectors, 
not  all  dy^^  are  achievable.  Therefore,  dy^*  may  need 
to  be  clipped  or  limited  such  that  u  =  p{dy^*)  does 
not  violate  actuator  limits.  Past  control  allocation 
research  includes  algorithms  that'  minimize  control 
deflections  [6]  or  drag  [9]  and  limit  unachievable 
commands  by  preserving  the  direction  of  the  orig¬ 
inal  command  [10,  11].  However,  it  is  shown  later 
that  preserving  the  direction  of  the  command  may 
undesirably  degrade  flying  qualities. 

The  following  control  allocation  algorithm  pro¬ 
vides  redundant  actuator  commands  that  optimize 
any  objective  that  is  a  function  of  the  control  deflec¬ 
tions.  The  algorithm  also  utilize  all  of  the  available 
control  power  and  limits  unachievable  commands  op)- 
timally. 


The  scalar  e  allows  the  algorithm  to  relax  the  equal¬ 
ity  constraint  if  dy®'  is  not  achievable  and  therefore 
provides  indication  whether  dy^*  is  achievable.  If 
6=1  then  dy^*  is  achievable,  otherwise  it  is  not. 
In  effect  c  limits  dy^\  if  necessary,  by  reducing  its 
magnitude  and  maintaining  its  direction.  So  this  op¬ 
timization  formulation  could  give  the  same  results  as 
some  past  control  allocation  work  [10,  11].  A  more 
general  control  allocation  framework  is  desirable  to 
provide  greater  flexibility  in  command  limiting.  Pri¬ 
oritized  control  leads  to  this  general  framework. 

3.3.2  Prioritized  Control 
Command  Limiting 

Recall  the  control  law  is  in  terms  of  a  desired  reduced 
dimension  control  dy^*.  Consider  decomposing  the 
control  law  into  a  set  of  partitions 

d**  =dl+<f2  +  .--  +  <ffc  (16) 

where  it  is  assumed  desirable  to  limit  dp‘,  if  nec- 
essary,  by  first  limiting  d^,  then  d^^i,  and  so  on. 
This  suggests  a  priority  among  the  control  law  par¬ 
titions.  Note  the  partitions  and  the  priority  among 
partitions  may  not  be  easy  to  determine  for  all  con¬ 
trol  laws.  A  possible  priority  for  a  dynamic  inver¬ 
sion  control  law  is  developed  in  a  later  section,  but 
for  now  assume  the  partitions  and  the  priority  are 
given.  The  following  constrained  optimization  prob¬ 
lem  ensures  the  priority 


3.3.1  Optimal  Objective  Control  Allocation 

Limits  on  the  controls  complicate  the  analytical  for¬ 
mulation  of  the  control  allocation  problem.  The  un¬ 
limited  control  zJlocation  problem  can  alternatively 
be  cast  into  the  following  linear  algebra  problem. 
Find  a  control  ti  such  that  Hyti  =  d^**  for  a  given 
dy*'.  To  exploit  redundemt  effector  suites,  the  fol¬ 
lowing  constrained  parameter  optimization  problem 
is  cast  from  the  linear  algebra  formulation 

nyn  J  =  /(«)  subject  to  ByU  =  d^®*  (14) 

Further  constraints  may  be  added  to  account  for  lim¬ 
its  on  u 


my  J  =  /(«)- 

subject  to 

jByii  —  Aidj  "b  ^2^2  "b  '  *  *  “b  Ajtdjt 
<  «<  <  ««i< 


0  <  A,-  <  1 


If  Aj  =  1  :  A,_i  =  Aj.2  =  •  •  •  =  1 
If  Aj  =  0  :  Aj'^i  =  Aj+2  =  •  •  •  =  0 

If0<  Aj  <  1 


.  /  Aj-i  —  Aj_2  —  • 

I  =  ^;+2  =  ■ 


I  =  1,  • ,  m 


(17) 


=  1 
=  0 


IQip  J  =  /(«)-€ 


subject  to 


( 


ByU  =  £<« 
0<€<1 


: «  =  1, 


(15) 


m . 


Prioritized  Dynamic  Inversion  Consider  the 
dynamic  inversion  control  law  with  proportional- 
plus-integral  compensation  in  eqs.(9)  and  (10).  The 
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3.3.3  Axes  Prioritization 


structure  of  this  control  law  lends  itself  to  the  fol¬ 
lowing  control  partitions 

=  da  H“  dp  +  d,*  +  dc  (18) 

da  =  — ty  — '  Ayyy  “I"  A^y 

Ad  =  diag(Ayy,,  ,  Ayy^^,  ’  *  *  ,  Ayy^^^^) 

dp  =  -Aay  -  wcy 
di  =  tJc^i 
dc  H  We/.y*"" 

where  Ayy,,.  is  the  (;,  j)  element  of  Ayy.  The  com¬ 
mand  part  (dc)  provides  aggressive  command  re¬ 
sponse  augmentation.  The  integral  part  (d,)  pro¬ 
vides  desired  loop  robustness  and  some  command 
response  augmentation.  The  proportional  part  (dp) 
provides  stability  augmentation.  The  cascading  and 
decoupling  part  (da)  provides  stability  and  decou¬ 
pled  response.  If  aggressive  reference  commands  are 
the  main  cause  of  unachievable  dy®^,  then  it  seems 
reasonable  to  limit  dc  before  d,-.  Since  basic  sta¬ 
bility  is  accomplished  by  da,  then  da  should  be  the 
last  part  that  is  limited.  Stability  augmentation  is 
accomplished  by  da  and  dp,  while  command  augmen¬ 
tation  is  accomplished  by  dc  and  d*.  Since  stability 
augmentation  is  most  critical,  the  following  priority 
is  proposed 

dy^*  =  Alda  +  A2dp  +  Aad,*  -f  A4dc  .  (19) 

Whereas  some  command  limiting  approaches  only 
reduce  feedforward  gains,  the  proposed  control  pri¬ 
ority  also  reduces  loop  gains  so  that  pilot  commands 
still  affect  the  aircraft  response. 

The  command  and  loop  responses  are  analyzed 
to  verify  the  proposed  priority.  In  each  output  chan- 


nel,  the  command  response  is  given  by 

y  XiUefeS  +  AaW^/i 

~  +  UeS  +  Xswlfi 

(20) 

and  the  loop  response  is 

s 

(21) 

It  is  seen  in  Fig.  (2)  that  the  loop  response  grace¬ 
fully  degrades  as  A3  goes  from  one  to  zero  followed 
by  further  graceful  degradation  as  A2  goes  from  one 
to  zero.  Similarly,  the  command  response  grace- 

fully  degrades  in  Fig.(3)  as  A4  goes  from  one  to  zero 
followed  by  further  graceful  degradation  as  A3  goes 
from  one  to  zero. 


Since  dy^*  is  a  vector,  further  priority  may  be  en¬ 
forced  between  its  elements  which  provides  axes  pri¬ 
oritization  capability.  For  brevity,  this  feature  will 
be  demonstrated  without  partitioning  dy®'.  How¬ 
ever,  axes  prioritization  capability  may  be  applied 
to  each  dy®'  partition  to  provide  the  most  general 
solution.  Assume  without  loss  of  generality  that  dy®' 
has  three  elements,  and  it  is  desired  to  limit  the  sec¬ 
ond  element  first,  the  third  element  second  and  the 
first  element  last.  The  following  scaling  matrix  will 
limit  dy*',  if  necessary,  according  to  the  stated  axes 
priority 

A  =  diag(Ai,  A3,  A2) .  (22) 


The  following  parameter  optimization  problem  is 
now  formulated  to  account  for  axes  prioritization 


;=i 

subject  to 
Byu  =  Ad** 

<  «•  <  «Ui 

0  <  A,  <  1 


If  Aj  =  1  :  Aj_i  =  Aj_2  =  •••  =  ! 
If  Xj  =  0  :  Aj+i  =  Aj+2  =  •  •  •  =  0 
>-i  =  =  • 

Aj+2  =  ■ 


IfO<  Aj  <  1: 


\ 


(23) 


To  combine  prioritized  control  with  axes  prioritiza¬ 
tion,  a  diagonal  matrix  (A^)  must  be  constructed  for 
each  control  partition  (dt)  to  reflect  the  prioritiza¬ 
tion  between  elements  of  that  partition. 


3.3.4  Implementation  Considerations 

Issues  primarily  related  to  computational  limits  be¬ 
come  critical  for  on-line  implementation  of  control 
allocation  algorithms.  If  the  control  law  is  to  be 
digitally  implemented,  control  limits  may  include  ac¬ 
tuator  position  limits,  rate  limits  multiplied  by  the 
digital  flight  control  system  sample  time,  and  any 
other  limits  that  may  be  discretized.  The  values  in 
the  optimization  problem  of  eq.(17)  should  become 
increments  from  current  values 

^  =  txo  +  «A  (24) 

dy^*  =  do  -F  dA 
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where  (  )a  indicates  increments  and  ( )o  current  val¬ 
ues.  For  sufBciently  fast  control  update  rates,  the 
optimization  objective  function  may  be  reasonably 
approximated  by  a  linear  function 

/(wo  +  «a)«C’J’(uo  +  «a)  (25) 

which  may  reduce  computational  requirements  of 
the  control  allocation  optimization.  However, 
the  linear  approximation  forces  boundary  solutions 
which  may  result  in  chattering. 


Control  Preference  Knowledge  of  a  control  pref¬ 
erence  and  adaptive  linear  objective  coefficients  may 
assist  in  the  elimination  of  chattering.  Assume  a 
control  preference  (u^)  that  does  not  violate  actua¬ 
tor  deflection  limits,  such  as  minimum  drag  or  mini¬ 
mum  wing  root  bending  moment  deflection,  is  given. 
In  general,  the  control  preference  does  not  achieve 
the  desired  y-derivative,  i.e.  ByUp  ^  The  so¬ 
lution  to  the  following  optimizatioh  problem 

nyn  J  =  !(«  -  UpfWf{u  -  Up)  (26) 
subject  to  5yti  =  do 

always  achieves  the  current  y-derivative  do  and  is 
given  by 

=  «p  +  W-^BE {ByW-^BFr\do  -  ByUf)  . 

(27) 

The  control  Uopt  is  “near”  tip  in  the  sense  of  the  op¬ 
timization  problem  in  eq.(26),  and  the  weight  Wp 
may  be  adjusted  to  guarantee  that  Uopt  will  not  vi¬ 
olate  actuator  deflection  limits.  The  linear  objec¬ 
tive  coefficients  in  eq.(25)  are  chosen  to  always  drive 
the  controls  toward  u^pi.  Once  a  control  effector 
is  within  a  tolerance  of  its  optimum  defined  by  the 
rate  limit  and  control  sample  rate,  it  may  effectively 
be  removed  from  the  optimization  by  adjusting  its 
coefficient  to  zero 


{ 


0  :  -  UoJ  <  tlrjAf 

-(liopt.  -  i/o<)  :  else 


where  tip  is  a  vector  of  rate  limits. 


(28) 


Sequential  Optimization  Since  the  constraints 
for  A  are  piece- wise  linear  in  eq.(17),  a  sequential 
approach  to  the  optimization  may  reduce  computa¬ 
tion  zJ  requirements.  The  idea  is  that  a  sequence  of 
k  optimizations  with  linear  constraints  and  a  sin¬ 
gle  A  element  is  solved  instead  of  a  single  optimiza¬ 
tion  with  nonlinear  constraints  and  k  elements  in  A. 


The  constraints  are  linear  for  each  sequential  opti¬ 
mization,  and  the  number  of  variables  is  reduced  by 
ib  —  1.  The  following  optimization  problem  is  solved 
initially  for  g  =  i 

^  =  /(«a)-A  (29) 

subject  to 

ByUt^  = 

<  «u, 

0<A<1 

1  =  !,•••, m. 

If  there  is  no  feasible  solution,  solve  with  g  =  t  —  1 
and  so  on.  The  solution  =  0,  A  =  0  is  always  fea¬ 
sible  since  the  current  control  satisfies  do  =  ByU^. 
So,  there  is  always  a  feasible  solution  in  this  sequence 
if  control  increments  (ti^)  are  used. 

Integrator  Windup  Dynamic  control  partitions 
may  experience  integrator  windup  for  coimnands 
that  are  unachievable.  However,  the  command  lim¬ 
iting  scalings  (A)  may  be  used  to  prevent  integrator 
windup.  For  example,  the  dynamic  inversion  inte¬ 
gral  control  partition  (d^)  in  eq.(19)  will  continue  to 
grow  due  to  the  Xi  dynamics  in  eq.(lO).  Integrator 
windup  is  prevented  by  stopping  the  integral  action 
and  holding  x,*  constant  whUe  A3  <  1. 

4  Tcdlless  Fighter 
Control  Allocation 

This  section  describes  application  of  optimal  objec¬ 
tive  control  allocation  to  a  tailless  fighter  aircraft. 
The  tailless  aircraft  is  a  65  deg.  sweep  delta  wing, 
single  engine,  mutli-role  supersonic  fighter  with  in¬ 
ternal  weapons  carriage.  A  high-fidelity  six  degrees- 
of-freedom  nonlinear  simulation  was  developed  from 
aerodynamic  wind-tunnel  data  of  a  scale  model  gen¬ 
erated  under  the  Innovative  Control  Effectors  (ICE) 
program  [3,  1].  The  ICE  configuration  includes  a 
large  suite  of  conventional  and  innovative  control  ef¬ 
fectors  that  provide  forces  and  moments  in  multi¬ 
ple  axes.  The  conventional  effectors  include  elevons, 
pitch  flaps,  thrust  vectoring,  and  outboard  leading 
edge  flaps.  The  innovative  or  unconventional  con¬ 
trol  effectors  include  spoiler-slot  deflectors  and  all- 
moving  tips. 
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A  dynamic  inversion  control  system  for  subsonic 
flight  has  been  developed  previously  for  an  ICE  con¬ 
figuration  without  flaps  and  spoiler-slot  deflectors 
[3].  The  past  work  includes  simple  ganging  control 
allocation,  and  the  current  work  develops  a  more 


sophisticated  control  allocation  method  for  the  ICE 
configuration  with  the  full  control  effector  suite. 

The  following  linear  model  is  generated  at  Mach 
0.4  and  15,000  ft.  altitude  for  demonstration  of  the 
optimal  objective  control  allocation 


or 

P 

p 

Q 

.  R  . 

L 

-0.6344 

0 

0 

-0.5971 

0 


0.0027 

0 

0.9871 

0 

-0.0038 

0.1540 

0 

-0.9876 

-8.2125 

-0.7849 

0 

0.1171 

0 

0 

-0.5099 

0 

-0.8887 

-0.0299 

0 

-0.0156 

-I- 


■  -0.0459 

-0.0459 

-0.0395 

-0.0047 

0.0047 

0 

3.7830 

-3.7830 

0 

-2.5114 

-2.5115 

-1.9042 

0.0453 

-0.0453 

0 

■  -0.0109 

0 

0.0217 

0 

0.0110 

0.0066 

0 

0.0790 

-2.0956 

-1.1329 

0 

1.5046 

0 

-0.8038 

-0.0283 

-0.0133 

-0.0133  ■ 

0.0031 

-0.0031 

1.8255 

-1.8255 

6pflap 

-0.9494 

-0.9494 

^amit 

-0.2081 

0.2081 

^amtr 

0.0217 

0.0047 

0.0047 

-0.0066 

-0.0021 

0.0021 

2.0957 

-0.3067 

0.3067 

1.5046 

-0.0003 

0.0004 

0.0283 

0.0937 

— 

0.0937 

-| 

Sptv 

Sytv 

^Stdl 

^$$dr 

Kuji 

(30) 


Note  that  this  model  is  of  the  form  in  eq.(8) 
without  the  bias  since  it  is  an  equilibrium  or  trim 
condition.  The  control  law  in  eqs.(9)  and  (10)  is 
completely  specified  by  the  uncommanded  variables 
(z),  commanded  variables  (y),  and  d^ign  parame¬ 
ters  (wc>  /ti  and  /c).  The  commanded  variables  are 
chosen  based  on  flying  qualities  requirements  [12]  to 
be  stability  axis  roll  rate,  pitch  rate,  and  a  blend 
of  sideslip  and  stability  axis  yaw  rate.  The  uncom¬ 
manded  variables  are  defined  to  be  angle-of-attack 
and  sideslip  which  guarantee  existence  of  the  trans¬ 
formation  inverse  in  eq.(4) 

[;] 

PcosQfo  +  ifeinoro 

Q 

—  Psinoo  +  iJcosao  _ 

-1. 

The  following  feedback  and  feedforward  gains  are 
chosen  based  upon  uncertainty  bounds  and  flying 
qualities  requirements  for  similar  fighter  aircraft  [6] 


sec 


(32) 


fi  =  25 

/e  =  .5. 

The  linear  model  of  the  aircraft  is  considered 
with  actuator  position  and  rate  limits  to  demon¬ 
strate  the  optimal  objective  control  allocation  and 
prioritized  control  concepts.  The  control  allocation 
optimization  objective  is  to  force  the  controls  to¬ 
ward  a  zero  control  preference  =  0),  i.e.  mini¬ 
mum  control  deflection.  This  is  accomplished  using 
control  increments  and  linear  objective  coeflBcients 
shown  in  section  3.3.4.  Although  the  algorithm  is 
capable  of  more  complex  objective  functions,  com¬ 
putational  requirements  may  increase  significantly 
for  complex  objectives.  This  may  be  important  for 
the  ultimate  goal  of  on-line  implementation. 

The  sequential  linear  program  in  eq.(29)  is  solved 
for  the  linear  tailless  fighter  aircraft  model  to  pro¬ 
vide  minimum  drag  control  allocation  for  achievable 
commands  and  prioritized  command  limiting  for  un¬ 
achievable  commands.  Note  that  there  are  eleven 
control  surfaces  so  m  =  11.  The  dynamic  inversion 
control  priority  is  indicated  by  the  following  con¬ 
troller  partitions 

di  =  drf  ^2  =  dp  da  =  di  =  dc  (33) 
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The  upper  and  lower  control  limits  are  defined  by  the 
most  restrictive  of  position  limit  (/p,  tip)  or  product 
of  rate  limit  (ur)  and  time  step  (Ai) 

ti/,  =  max{/p,- -  tio,, -tirjA<}  (34) 

Uu,  =  min  {tip,- —  tiOij  tirj At} 

where  uq,  Is  the  element  of  the  current  control 
uq.  The  sequential  linear  program  is  solved  at  each 
simulation  time  step  for  the  following  analysis, 

4.1  Simulation  Analysis 

Simulation  responses  to  simultaneous  step  com¬ 
mands  in  roll  (yi)  and  pitch  (y2)  are  presented  in 
this  section  to  show  the  benefits  of  the  prioritized 
control  command  limiting.  A  step  pitch  command 
at  1  sec.  is  followed  by  a  step  roll  command  at  3  sec. 
for  two  cases.  The  first  case,  indicated  in  Figs. (4) 
and  (5)  by  prior,  is  optimal  objective  control  alloca¬ 
tion  with  prioritized  control  command  limiting  (see 
eq.(17)).  The  second  case,  indicated  in  Figs. (4)  and 
(5)  by  dir,  is  optimal  objective  control  allocation 
with  command  limiting  that  maintains  command 
direction  (see  eq.(15)).  The  command  variable  re¬ 
sponses  are  shown  in  Fig. (4),  and  the  scalings  for 
unachievable  commands  are  shown  in  Fig.  (5).  Note 
that  dy*'  is  unachievable  at  the  onset  of  the  com¬ 
mands  as  indicated  by  e  <  1  and  A4  <  1.  Fur¬ 
ther,  note  the  coupling  of  the  roll  response  with  the 
pitch  and  yaw  response  due  to  limiting  commands  to 
maintain  command  direction.  However,  the  priori¬ 
tized  control  command  limiting  only  limits  dc  which 
maintains  the  desired  decoupled  responses.  Also,  the 
dir  case  rise  time  is  longer  for  the  roll  response.  The 
control  responses  are  available  but  are  not  shown 
due  to  space  limitations. 

4.2  Future  Plans 

Loop  stability  during  command  limiting  may  be 
guaranteed  for  certain  fixed  scalings  (A),  but  the 
scalings  are  generally  dynamic.  Stability  analysis 
that  captures  the  dynamic  nature  of  the  scalings  of 
the  prioritized  control  command  limiting  is  planned. 
To  date,  the  current  algorithms  have  only  been 
demonstrated  on  linear  airframe  models  with  nonlin¬ 
ear  actuators.  Future  plans  include  extension  of  op¬ 
timal  objective  control  allocation  to  the  ICE  tailless 
fighter  nonlinear  6DOF  simulation.  The  nonlinear 
6D0F  will  allow  further  analysis  such  as  computa¬ 
tional  requirements,  mode  transition  effects,  robust¬ 


ness  to  unmodelled  dynamics  and  parametric  uncer¬ 
tainty,  and  effects  of  additional  control  effector  con¬ 
straints.  Efforts  to  minimize  computational  require¬ 
ments  of  the  algorithms  have  been  accomplished,  but 
a  comprehensive  computational  requirements  analy¬ 
sis  is  necessary.  Analysis  of  transients  associated 
with  control  allocation  mode  transition,  from  mini¬ 
mum  drag  to  minmum  wing  load  for  example,  will 
be  evaluated.  Some  noise  and  uncertainty  has  been 
injected  into  the  simulations,  however  thorough  ro¬ 
bustness  analysis  is  planned.  Control  effector  inter¬ 
actions,  hinge  moment  limit  effects,  and  loads  limit 
effects  will  be  analyzed  in  the  future  as  well. 

Other  plans  include  using  a  modified  sequential 
least  squares  identification  algorithm  [13, 14]  to  pro¬ 
vide  stability  and  control  derivatives  to  the  augmen¬ 
tation  and  control  allocation  modules  in  real-time. 
This  will  give  the  control  system  reconfiguration  ca¬ 
pability  in  the  event  of  battle  damage  or  component 
failure.  The  reconfiguration  aspect  of  tailless  fighter 
flight  control  is  investigated  in  the  RESTORE  pro¬ 
gram  [Ij.  Control  redundancy  will  no  doubt  compli¬ 
cate  real-time  control  derivative  identification.  How¬ 
ever,  the  optimal  objective  control  allocation  algo¬ 
rithm  may  possibly  be  augmented  to  assist  in  the 
identification  of  redundant  effector  control  deriva¬ 
tives. 

5  Conclusions 

Optimal  objective  control  allocation  with  prioritized 
control  command  limiting  is  developed  as  a  general 
framework  for  redundant  effector  management  and 
optimal  command  limiting  to  prevent  actuator  sat¬ 
uration.  The  control  allocation  problem  is  posed 
as  a  constrained  optimization  problem  with  added 
variables  for  command  limiting.  Although  the  al¬ 
gorithm  is  valid  for  complex  objectives,  linear  ob¬ 
jectives  lead  to  linear  programs  which  may  reduce 
computational  requirements.  The  control  allocation 
algorithm  has  been  successfully  demonstrated  on  a 
tailless  fighter  aircraft  to  minimize  a  linear  approxi¬ 
mation  of  drag.  Prioritized  dynamic  inversion  com¬ 
mand  limiting  maintains  decoupled  axes  responses 
for  unachievable  commands,  while  past  command 
limiting  strategies  show  unacceptable  axes  coupling. 
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Figure  1:  Flight  Controller  Structure 
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Figure  4:  Command  Variable  Responses 
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Figure  5:  Command  Scaling  Responses 
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